Energy dependence of current noise in superconducting/normal metal junctions 
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Interference of electronic waves undergoing Andreev reflection in diffusive conductors determines 
the energy profile of the conductance on the scale of the Thouless energy. A similar dependence 
exists in the current noise, but its behavior is known only in few limiting cases. We consider a 
metallic diffusive wire connected to a superconducting reservoir through an interface characterized 
by an arbitrary distribution of channel transparencies. Within the quasiclassical theory for current 
fluctuations we provide a general expression for the energy dependence of the current noise. 

PACS numbers: 74.45.+c,74.40.+k,72.70.+m 



Interference of electronic waves in metallic disordered 
conductors is responsible for weak localization correc- 
tions to the conductance If these are neglected, the 
probability of transferring an electron through the diffu- 
sive medium is given by the sum of the modulus squared 
of the quantum probability amplitudes for crossing the 
sample along all possible paths. This probability is de- 
noted as semiclassical, since quantum mechanics is nec- 
essary only for establishing the probability for following 
each path independently of the phases of the quantum 
amplitudes. In superconducting/normal metal hybrid 
structures, interference contributions are not corrections, 
they may actually dominate the above defined semiclassi- 
cal result for temperatures and voltages smaller than the 
superconducting gap. This is seen experimentally as an 
energy dependence of the conductance on the scale of the 
Thouless energy. Indeed, the energy dependence comes 
from the small wavevector mismatch, linear in the energy 
of the excitations, between the electron and the Andreev 
reflected hole. This is responsible for the phase difference 
in the amplitudes for two different paths leading to inter- 
ference. The effect is well known and explicit predictions 
and measurements exist for a number of systems 

Interference strongly affects the current noise too 0. 
The largest effects are expected in the tunneling limit, 
when the transparency of the barrier is small and its re- 
sistance is much larger than the resistance of the diffusive 
normal region. Then, the conductance has a strong non 
linear dependence at low bias (reflectionless tunneling) 
[HQ. This is actually the case, but the zero-temperature 
noise (or shot noise) does not give any additional infor- 
mation on the system since it is simply proportional to 
the current, as shown numerically in a specific example 
in Ref. |5| and quite generally in Ref. 0. The double 
tunnel barrier system has been considered in Ref. || . In 
the case of a diffusive metal wire in contact with a super- 
conductor through an interface of conductance Gb much 
larger than the wire conductance Gd, Belzig and Nazarov 
9] found that the differential shot noise, dS/dV, shows 



a reentrant behavior, as a function of the voltage bias, 
similar, but not identical, to the conductance one. (The 
extension of the Boltzman-Langevin approach to the co- 
herent regime in Ref. neglects this difference.) In 
order to compare quantitatively with actual experiments 
[ill flil fl^ and to gain more insight in the interference 
phenomenon, it is necessary to obtain the energy depen- 
dence of noise in moregeneral situations. The numerical 
method used in Ref. is, in principle, suitable to treat 
more general cases, notably the case when Gd > Gb-, but 
only if all channel transparencies, {r„}, that characterize 
the interface are small. When arbitrary transparencies 
are present one has to solve numerically an additional 
self-consistent equation 01 . This appears particularly 
heavy numerically if one is interested to treat a distri- 
bution of transparencies. We are not aware of results in 
this direction. In this paper we present an analytical so- 
lution for the diffusion-type differential equation for the 
noise within the theory of current fluctuations |l5j in the 
quasiclassical dirty limit HJ. It allows to treat the gen- 
eral case of arbitrary values for {r„} and Gb/Gd- We 
express the noise in terms of a function that satisfies a 
linear differential equation to be solved numerically once 
the channel distribution is given. We can thus isolate the 
energy-independent (semiclassical) contribution and the 
interference contribution to both the conductance and 
noise. We discuss their relation in the following. 

We begin by stating the framework for the theory 
of current fluctuations 0, ^J. It relies on the eval- 
uation of the quasiclassical Green's functions g(x,e,x) 
in the Nambu(")-Keldysh(~) space, at a given energy e 
and counting field \. The counting field appears as a 
gauge transformation of the Green's function in the nor- 
mal reservoir 



9n(x) 



(1) 



where g% = f 3 ® a 3 + (/to + fwh) ® (<ti + icr 2 ) is 
the standard Green's function for a metallic reservoir, 
o^i,Tj^j =1 2t 3) are Pauli matrices, f K = f 3 <E> a x , f L0 = 



2 



1 -/+-/_, /to = /- - /+, f±(e) = f(e ± eV), / is the 
Fermi function at temperature T, and V is the voltage 
bias. Given g(x), one can calculate the spectral matrix 
current in the wire 



J(x) = -LG D g{x)d x g{x) 



(2) 



with x the coordinate along the wire of length L. The 
current and the zero-frequency noise are given respec- 
tively by @ 



I = J(X = 0) and S 



-He 



dJ(x) 



dx 



(3) 



where J(%) = — l/(8e) Jde Tr [fx J(x)] and J(x) does 
not depend on the position x where J is evaluated. 

The Green's function g(x) in the wire is determined by 
the diffusion-like Usadel equation 



HD 



d x J(x) = -is [f 3 ,g(x)] 



(4) 



(where f 3 = f 3 <g> 1 and D is the diffusion constant in the 
wire) and by boundary conditions at the two extremities. 
We assume a good contact on the normal side at x = 
L. This implies g(L) = g^ as defined in On the 
superconducting side, at x = 0, the boundary condition 
expresses the conservation of the spectral matrix current 
© through the interface [3,E3 : 



2T n [g(0),g s } 



4 + r„({5(0),.g s }-2) 



(5) 



where the eigenvalues T n of the transmission matrix 
through the interface appear explicitly and Gb = 
Gq Y] r„, with Gq = e 2 / (ttH) the quantum of conduc- 
tance. The Thouless energy, Et = HD/L 2 , is the relevant 
energy scale in Eq. Q. We restrict to the case where it 
is much smaller than the superconducting gap. With the 
same restrictions on the voltage bias and temperature, 
the Green's function in the superconducting reservoir is 
9S = T2 <8> 1. Extension to higher energies is straightfor- 
ward, but it will not be considered to keep the presen- 
tation simple. The quasiclassical Green's functions obey 
the normalization condition g 2 = 1 and the symmetry 
property 



g(x)^ = tl g(x) tl, with 



TL = T 3 (g) (72 



(6) 



Finding J(x) f° r an values of x is equivalent to calcu- 
late the full counting statistics of charge transfer 0| and 
it may be a formidable task. One of the main difficulties 
comes from the normalization condition g = 1. As a 
matter of fact, for x — the Green's functions have a 
triangular structure in Kcldysh space: 



( 9 Q f A 



(7) 



In the absence of supcrcurrent it exists a simple 
paramctrization fulfilling the normalization condition 



and the symmetry property 



9 



if 2 sinh6», g K = g R f - fg A , and f = ft 



f 3 cosh 9 + 
Tl -I- r 3 f T , with 
9 = 6! + i9 2 and f L , f T , 9 X , and 9 2 real. When X + 
the triangular structure is lost and we are not aware of a 
simple paramctrization for g in this case. 

A less ambitious program is to calculate only the noise. 
In this case it is clear from Eq. that the full depen- 
dence of g(x) is not needed: dg(x)/dx\ x =o is enough. As 
suggested in Ref. [t| one can thus develop g(x) i n powers 
of X- 



g(x, x) = go(x) - i(x/2)gi(x) + 0( x 2 



(8) 



and solve the problem order by order in x- F° r the lead- 
ing order Green's function, go, the above paramctrization 
is appropriate. Eq. (@J for g R gives the non-linear equa- 
tion 



hD9"(x) + 2iesmh 9 (x) = 0. 
with boundary conditions 9(L) = and 



L9'(0) 



cosh (9(0) 



f 



(isinh0(O) - 1) 



(9) 



(10) 



We defined (^(T)) = £ n r«V>(T n )/ £„ r « and r = 
Gd/Gb- Eq. 0} for g K gives Jl{x) = fw and the dif- 
ferential equation for Jt'- (cosh 2 9i(x)f T (x))' = with 
boundary conditions Jt(L) = fxo an d 



/t(0) 



St (O)0i(O) 



cosh 6»i (0) sinh0i(O) 
Then the current is / = l/(2e) J(kQ(e)f T o{e) with Q 



(11) 



G(e) = G L 



tanh6»i(0) 
L9[(0) 



(12) 



and V~ 1 (e) = 1/L j Q L ds/ cosh 2 #i(s). At low tempera- 
tures, k B T < eV, G(V) = G(eV) is the differential con- 
ductance. These equations have been recently exploited 
in Ref. 01 to discuss the conductance. 

At the next order in x, we get a linear matrix differ- 
ential equation for g\{x): 



HDd x J 1 (x) =ieLG D [f 3 ,9i{x)] 



(13) 



with Ji(x) = -LG D [go{x)d x g 1 (x) + gi(x)d x g (x)}. The 
boundary conditions read gi(L) = on the nor- 

mal side and Ji(0) = 2G B {ABA) where A = [4 + 
r({ff(0),ffs} - 2)]- 1 and B = [4 + 2T(g s9 o(0)gi(0)gs ~ 
go{0)gi{0) — [51(0), gs})] on the superconducting side. 
Finally, the normalization of g gives the condition 
{go{x), gi(x)} = that can be fulfilled with the change 
of variable gi{x) = [go(x), <fi(x)]. The matrix <j)(x) is now 



3 



constrained by the symmetry properties JBJ only. We 
find that Eq. (|13f) can be conveniently solved with the 
following parametrization for </>: 



a /rem - c / f 3 bf 3 + d 

c r 3 a* / T0 fi + c / f 3 



(14) 



with a = ai + ia2 and ai, a2, 6, c, d real functions of 
x. Substituting g\ in terms of 4> into (|13|l after straight- 
forward but lengthly calculations we obtain the set of 
equations and boundary conditions for the four parame- 
ters a, b, c, and d. 



The parameter a plays the role of in the lowest order equation 10 : 

sinh0i(x) G{e) 2 



HDa"(x) +2iea(x) cosh0(x) = -2E T 



cosh 3 0i (x) G\ 



(15) 



The boundary conditions are a(L) = and L a'(0) = a a(0)/r + /3/r with 



Q 



isinh0 - T(isinh0 - l)/2 
[l + r(isinh0 - l)/2] 2 



/3 



= 2 / 2r 2 cosh 0* + 8 (r - 1 ) cosh - 2ir(r - 2) sinh cosh 0* 



|l + r(isinh0 - l)/2r (1 
I 



r(isinh0- l)/2) 



(16) 



both evaluated at x = 0. The equation for b resembles 
that for fx- it is solved analytically in terms of an inte- 
gral that enters the final expression for the noise. The 
parameter c turns out to be simply proportional to fr(x): 
c(x) = -/t(x)// to with c(0) = 1 - g(e)/[G D V(e)]. Fi- 
nally we find d(x) = 2/ lo /to[1 + c 3 (x) — 2 tan 02 (x)]. Its 
knowledge is not necessary to obtain the noise. 

The expression for the noise is obtained by evaluating 
Eq. ©: 



deG(e) {1 - flo(e) [1 - F(e)]fh(e)} , (17) 



where 



He) = g(l 



-c(0) 3 ) 



2G(e 



G 



D 



,(0) 



Gz,oi(0)c(0) 



Z" 1 sinh0iai 
/o cosh 3 0i 
2oi(0) 



5(e)tanh0i(O) sinh20i(O) 



(18) 



In equilibrium, Eq. (|17|l yields the fluctuation dissipation 
relation S = iksTG, as expected. Out of equilibrium, 
in the shot-noise regime ks T eV, Eq. (|18|) defines the 
experimentally accessible differential Fano factor F(V) = 
(dS{V)/dV)/(2eG(V)) =F(eV). 

Eq. (|18|l is the central result of the paper. Once the 
problem for the conductance has been solved and 0(x) 
and /t(x) are known 0,^3' we provide a simple and ef- 
ficient way to calculate the noise. It suffices to solve the 
linear differential equation (|15fl with the given boundary 
conditions, and substitute the result into Eq. (|18|) . This 
program has to be followed numerically in most situa- 
tions, but its implementation is straightforward and al- 
lows to obtain quantitative predictions for a wide range 
of realizable experiments. The previous approach to the 
same problem has been to discretize Eq. (0J from the 
outset and solve it numerically 0, 0] . That technique 



is particularly efficient to treat the simplified boundary 
conditions J(0) = 2Gb [<?(0), gs] E3i valid only when all 
r„ <C 1. Its application to the general case is numerically 
much more difficult. 

We can now discuss the crossover from the coherent 
(e <C Et) to the incoherent semiclassical (e 3> Et) 
regime quite generally. 

In the fully coherent regime we recover known expres- 
sions obtained with random-matrix theory [5ll20j| or qua- 
siclassical Green's functions [U- In the opposite limit, 
propagation is incoherent: due to the large value of e the 
phases accumulated along two different paths are always 
uncorrelated. In order to verify this point we first re- 
peated the calculation when both reservoirs are normal. 
It suffices to substitute gs with gx into (JSJ with eV = 0. 
(The procedure greatly simplifies since g becomes diago- 
nal in Nambu indexes.) As expected, we find no energy 
dependence for G and F. For the conductance the usual 



Ohm's law holds: G 
factor we find 



G 



D 



G R while for the Fano 



F 



1 



2-3 



y r 2 
E„r„; {g d + g b y 



r 3 



(19) 



Eq. (|19fl coincides with the semiclassical result obtained 
in Rcf. . In the superconducting case, Eq. i|12|) 
and Eq. i|18[) can be evaluated analytically for e ^> Et 
defining G nlc and F; nc , respectively. It turns out that 
the conductance and the Fano factor are given by the 
same expressions for the normal case multiplied by a 
factor 2 and where the substitutions Go — * Gd/2 and 
r„ — > = T 2 /(2 - r n ) 2 have been operated (the last 
one implies Gr -> Gr <r/(2 - T) 2 )). This result can be 
understood jU I24I l25| in terms of semiclassical proba- 
bility for charge transfer, that is the probability for an 
electron to cross the distance L to the barrier, to be An- 
dreev reflected, and then to come back as a hole on the 
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FIG. 1: Energy dependence of the conductance and of the 
differential Fano factor when the wire and the interface have 
the same conductance (Gd = Gs) for different values of V = 
F n . The dashed lines correspond to the universal distribution 
{T n } for disordered interfaces. 



same distance. Therefore, the wire in contact with the 
superconductor can be described by two normal wires, 
one for the electrons and the other for the holes, coupled 
through a fictitious barrier that converts electrons into 
holes with probability and to which a voltage bias of 
2 V is applied [13 . 

Quantum corrections to both G(e) and F(e) can be 
evaluated as an expansion for large s. We give the re- 
sult for_all T„ = T: G = G inc + Gi(e) , with Gi = 
2 y/E T /e (4 - 4T 



r 2 )/(r(r-2) 2 + 2r) 2 , and 
F inc -r-y(T,r)[G(e)-G inc \ 



(20) 



with 7 = [r(2 - T) 2 (64 + (2 - T) T (64 + r (2 - T) (12 - 
12T - T 2 )))]/[(r(2 - T) 2 + 2 T) 2 (A- AT- T 2 )]. 7 is posi- 
tive for most values of T. Eq. 1121 III relates the interference 
contributions to G with those to F for large e. When the 
barrier dominates (r ^ 1) the relation (|2L)|I holds for any 
value of e up to first order in 1/r with j(T) = "f(T, r — > 
00) and G^e) = 20(2^^7^) (4- 4r - T 2 )/[r 2 (T-2) 4 ] 
and (j)(x) = (sinh.T + sin 2;) / [2 (cosh 2; + cosx)]. This sug- 
gests that the simple relation (f^TTf) stands, though ap- 
proximately, beyond the range of validity for which it 
has been proved. A possible interpretation is that inter- 
ference modifies the effective transparency of the whole 
system. The linear relation between the quantum correc- 
tions to G and F then corresponds to the single channel 
quantum result: G oc T and F oc 1 — T 

For intermediate values of the parameters r, e, and T n , 
quantum interference contribution can be studied numer- 
ically. We report the results in Fig. 1 for Gd = Gb and 
different values of T = T n for all n. The transparency 
r drives a smooth crossover from a reflectionless tun- 
neling to a reentrant behavior which can be seen both 
in the conductance and the Fano factor. This proves 



that the energy dependence of both may be strongly af- 
fected by the precise set {r n } of transparencies of the 
interface. As a relevant example we plot the result 
for the universal distribution for disordered interfaces 
En 6 ( T ~ r «) « I/^a/I - T) M- Fig. 1 also con- 
firms that Eq. (|20|l is actually qualitatively satisfied in 
the whole range of parameters. 

In the limit Gd <C Gb (r — > 0) the actual values 
of {r„} drop from the boundary conditions. Indeed 
Eqs. tnnnHSl force 6>(0) = -in/2, f T (0) = (and 
thus c(0) = 0), and a(0) = 0. The conductance be- 
comes Q(e) = V{e)Gd and the Fano factor is given by 
Eq. I|18|) where only the 2/3 and the integral terms sur- 
vive. The first term is the semiclassical incoherent value 
[2?| . which also coincides with the fully coherent one [jj. 
The integral term singles out the interference contribu- 
tion to the Fano factor. The energy dependence coin- 
cides with that obtained previously in a different way 
which, in turn, agrees qualitatively with the experi- 
mental result |l2j . However, a broader voltage range for 
the reentrance in G(V) and S(V) is predicted j28|. A 
possible explanation is that the barrier resistance is not 
negligible. Assuming a disordered interface, we find that 
the reasonably small value of r = 0.3 allows to fit the 
conductance. It also improves the fit for the noise, but 
the agreement is not perfect. This may be due either to a 
different distribution of transparencies at the barrier, or, 
as suggested in Ref. 0, to heating or interaction effects. 

In conclusion, we provide a framework to calculate 
both conductance and noise in a normal metallic wire 
connected to a superconducting lead through an arbi- 
trary interface. We predict a strong energy dependence of 
the coherent contribution to the Fano factor. We suggest 
to exploit this dependence to experimentally characterize 
the transparency of interfaces. 
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